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Invariant circles play an important role as barriers to transport in the dynamics of area- 
preserving maps. KAM theory guarantees the persistence of some circles for near-integrable 
Ch maps, but far from the integrable case all circles can be destroyed. A standard method for 

1—t determining the existence or nonexistence of a circle, Greene's residue criterion, requires the 

I " I computation of long-period orbits, which can be difficult if the map has no reversing symmetry. 

We use de la Llave's quasi-Newton, Fourier-based scheme to numerically compute the conjugacy 
^~~^ of a Diophantine circle conjugate to rigid rotation, and the singularity of a norm of a derivative of 

i^ the conjugacy to predict criticality. We study near-critical conjugacies for families of rotational 

1— invariant circles in generalizations of Chirikov's standard map. 

l/-v A first goal is to obtain evidence to support the long-standing conjecture that when circles 

f~ — breakup they form cantori, as is known for twist maps by Aubry-Mather theory. The location 

_il of the largest gaps is compared to the maxima of the potential when anti-integrable theory 

f^ applies. A second goal is to support the conjecture that locally most robust circles have noble 

C^ rotation numbers, even when the map is not reversible. We show that relative robustness varies 

^"^ inversely with the discriminant for rotation numbers in quadratic algebraic fields. Finally, we 

J> observe that the rotation number of the globally most robust circle generically appears to be a 

piecewise-constant function in two-parameter families of maps. 



X 



1 Introduction 

Area-preserving maps are discrete-time analogs of Hamiltonian systems and represent the simplest 
of such systems to exhibit chaotic dynamics [Mei92]. A well-studied family are the generalized 



* Adam.Fox@colorado.edu 
'James.Meiss@colorado.edu 

'"AMF and JDM and were supported in part by NSF grant DMS-1211350. Useful conversations with Renato 
Calleja, Rafael de la Llave, Keith Julien, Robert Easton, and Holger DuUin are gratefully acknowledged. 



standard maps, f^-.M^ M, where M = 'E> xM., given by 

x' = X + Q-{y + £g{x)) mod 1, 

y = y + £9{x). 

We assume that the force, g{x), is a smooth, periodic function of the angle x with period one and zero 
average, so that fe has "zero net flux." The function Q is the frequency map, its dependence upon 
the momentum y represents the anharmonicity of the dynamics. When g and O are differentiable, 
(1) is an area-preserving diffeomorphism. We will assume that g and are analytic. 
Most studies of the dynamics of (1) use the force 

51(2;) = — sin(27r2;), (2) 

and the frequency map 

^i{y) = y- (3) 

We will refer to this choice as Chirikov 's standard map [Chi79] . This map applies to many physical 
models, including the dynamics of a kicked rotor, a charged particle in electrostatic waves, the 
adsorption of a layer of atoms on a crystal surface, and the motion of a particle in a relativistic 
cyclotron, see the references in [Mei92]. 

When £ = 0, all of the orbits of (1) lie on invariant circles, Ty = S x {y}, and the dynamics on 
each circle is simply rigid rotation with rotation number co = 0(yo), i.e., {xt,yt) = {xo + ujt, yo). We 
refer to any circle that is homotopic to To as rotational. Since the restriction of fs to a rotational 
invariant circle is a degree-one homeomorphism, each orbit on the circle has the same rotation 
number uj. 

KAM theory predicts that invariant circles with "sufficiently irrational" ^ rotation numbers per- 
sist when £ is sufficiently small [dLLOl, PosOl]. Alternatively when e is sufficiently large, converse 
KAM theory predicts that there are no rotational invariant circles [MP85]. The study of the 
breakup of these circles was initiated in the work of Greene [Gre79] and led to the development of 
a renormalization theory of breakup [Mac93] . 

The standard method for studying this breakup is Greene's residue criterion, recalled in App. B. 
Its implementation requires accurate computation of periodic orbits and the linearization of the 
map about these orbits. This computation is much easier when the force g is odd. In this case, 
(1) is reversible: there is an orientation-reversing involution that conjugates fe to its inverse." The 
fixed sets of this reversor intersect every rotational circle, see App. C. A consequence is that every 
rotational circle is symmetric and can be approximated by a sequence of symmetric periodic orbits. 
Symmetric orbits can be found by a one- dimensional root-finding algorithm. Due to this additional 
structure, all previous studies of breakup have assumed odd forcing. 

In this paper we will study the breakup of invariant circles for the generalized standard map 
when it lacks such symmetries. The mechanism for this breakup in nonreversible maps remains 
an open question, first posed by MacKay in his 1982 thesis, reprinted in [Mac93]. Our goal is to 
provide evidence for three conjectures. 

Conjecture 1. When a rotational invariant circle of (1) is destroyed it becomes a cantorus. 



'^ Namely, u) is Diophantine, see App. A. 

■^Since f2i is odd, Chirikov's map also has a second reversor, so we call it a "doubly reversible" map. 



This is known to be true for nondegenerate twist maps, i.e., when Drt{y) ^ 0. For this case, the 
powerful theory of Aubry and Mather iniphes the existence of recurrent "minimizing" invariant sets 
for each rotation number oj [Mat82, Aub83]. These invariant sets always lie on Lipschitz graphs 
over X, and when oj is irrational they are either rotational invariant circles or invariant Cantor sets, 
Percival's cantori [Per 79]. Thus a Diophantine invariant circle for a twist map persists up to a 
critical parameter value, Ecr, where it loses smoothness, and then for larger e becomes a cantorus.'^ 
The notion of "minimizing" arises from the Frenkel-Kontorova interpretation: orbits of a twist 
map are equivalent to equilibria of a one- dimensional, nearest-neighbor-coupled chain of particles. 
For example, define the periodic external potential V by DV{x) = g{x), and, for (3), the spring 
potential energy ^(x' — x)^ for each pair of neighboring particles at positions x, x' E M, to give the 
energy 

E{x,x') = l{x' -xf + eV{x). (4) 

Any sequence {xt : t G Z} that is a critical point of the (formal) sum £ = ^^£'(xi,Xi+i) is 
equivalent to an orbit of (1) for the frequency map (3). Aubry-Mather sets are minima of the 
energy with respect to variations with compact support. 

The creation of the gaps in a cantorus is most easily understood through the anti-integrable 
(AI) limit [AA90]. As e increases, the potential energy in (4) will begin to dominate the spring 
energy and particles in the minimizing state will tend to fall into potential wells. For large enough 
e, there will be no particles in a neighborhood of the maxima of V, opening gaps in the circle. The 
orbits of each gap form a bi-infinite family that Baesens and MacKay called a hole. More formally, 
AI theory studies the continuation of the critical points of lime_>oo -£, to finite e. In the limit, the 
particles simply sit on any sequence of critical points of V. If these critical points are nondegenerate 
and the "acceleration" of the sequence is bounded, then such states may be continued to finite, 
large enough e. Ordered AI states with irrational u continue to cantori [MM92]. 

However, neither Aubry-Mather theory nor the Frenkel-Kontorova energy (4) apply when the 
frequency map does not satisfy the twist condition. Thus, it is not known whether invariant circles 
of nontwist maps are destroyed by the formation of a family of gaps, nor whether there are remnant 
Cantor sets that remain. 

One well-studied example is the so-called standard nontwist map of Howard and Hohs [HH95], 
where the frequency map is 

n2{y,6) = y^-6. (5) 

KAM theory still applies to this map, providing — as for (5) — the curvature of Q does not van- 
ish when the twist vanishes [DdlLOO]. The nontwist map has rotational invariant circles that are 
meandering, in the sense that they are not graphs over x, and the nontrivial fixed point of the renor- 
malization group for critical circles in this map has been extensively studied [AWM03, AdlLPOS, 
FWAM06, WFPll]. 

In §4 we will study the formation of gaps in invariant circles of generalized standard maps with 
a variety of frequency maps and multi-harmonic forces, giving evidence for Conj. 1 in asymmetric 
and nontwist maps. 

Greene conjectured that the last invariant circle of Chirikov's map has the rotation number 
7 = ^(1 + \/5), the golden mean. Though this conjecture has never been proved, it is supported by 

^Of course, it may reform for larger values of e, and indeed the boundary of existence of an invariant circle for 
multi-parameter maps is often quite complex [Bul86, Wil87, KM89]. 



strong numerical evidence [Mac93]. More generally the golden invariant circle will not necessarily 
be globally most robust, but one expects that the last invariant circle in any neighborhood has a 
related rotation number: 

Conjecture 2. Rotational invariant circles with noble rotation numbers are locally most robust. 

A frequency ratio oj = 1^1/1^2 is noble if the vector z^ is an integral basis for the quadratic field 
Q(7), or equivalently if its continued fraction has a golden mean tail, see App. A. Strong numerical 
evidence for the local robustness of the nobles for Chirikov's map was given by MacKay and Stark 
[MS92]. In §5, we will study the relative robustness of circles with rotation numbers from different 
quadratic fields for several examples of (1). 

The local robustness conjecture implies that the critical function £cri^) will have local maxima 
at each noble rotation number. The global maximum of this function for Chirikov's map occurs, 
according to Greene, at the golden mean. More generally, let uj^ax be location of the global 
maximum of Sen i-e., the rotation number of the globally most robust circle. We will investigate 
the dependence of comax on the force g. 

Conjecture 3. The rotation number of the globally most robust invariant circle is generically 
piecewise constant under continuous changes of (1). 

We are not aware of previous studies of this conjecture. Our preliminary evidence consists of 
studying two-parameter families fe,^, obtained by introducing an additional parameter ?/) into g. 
The critical function will of course depend upon tp and if it were a typical, smooth function, the 
location of its global maximum would be a piecewise continuous function of V' with occasional 
jump discontinuities as different local maxima take over as the global maximum. As we will see 
in §6, ujmax{tp) seems instead to be locked to fixed noble rotation numbers over intervals of ■;/'• 
This conjecture does not hold when the map has special symmetries that cause ecr(w,'0) to be a 
discontinuous function of -0 for fixed cj. 

Each of these studies requires a technique for computing invariant circles of the generalized 
standard map that does not rely upon finding symmetric periodic orbits. We will use, see §2, a 
method introduced by de la Llave and collaborators [dlLGJVOS] to compute the Fourier series for 
the embedding k : S ^ M that conjugates the dynamics of fs on an invariant circle to the rigid 
rotation, T^^ : S — )• S, 

Tue) = e + u. (6) 

This technique gives a numerically robust method for computing subcritical invariant circles. We 
will then estimate the critical value by continuation in e and extrapolating to the critical point by 
estimating where the circle loses smoothness. In particular, a Sobolev norm of the conjugacy is 
singular at criticality [AdlLPOS, OP08, CdlLlOa, CdlLlOb]. We will see in §3 that a certain Sobolev 
seminorm has a power law blowup as e — )• £cr that allows accurate estimates of the critical value. 

2 Computing the Conjugacy for an Invariant Circle 

In this section we briefly recall an efficient algorithm to compute the embedding for invariant tori 
of maps, adapted from one developed by [dlLGJVOS, FdlLSOQ]. We specialize to the case that / 



is an area-preserving map and the torus is a rotational circle, though much of the analysis carries 
over to higher-dimensional symplectic maps. 

Suppose that / is a real analytic, area-preserving map of the cylinder M = S x M. We assume 
that / has an analytic, rotational invariant circle C C M on which the dynamics is conjugate to 
rigid rotation with Diophantine rotation number w, that is, we suppose there exists an analytic 
embedding A; : S — )• M such that k(S) = C, and 



f o k = koT^ 



(7) 



where T,^ is the rigid rotation (6), see Fig. 1. 
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Figure 1: A visualization of the conjugacy (7) and of automatic reducibility. Here the u is the unit tangent vector 
to the embedded circle C, v is the unit normal, and the shaded regions have unit area. 

Lifting the map / to the universal cover of M, the condition that C be rotational is 

k{0 + i) = k{e) + {i,of, 

i.e., the angle component of /c is a degree-one circle map (called the "hull" function in [Aub83]) 
and the action component is periodic. We will say, in this case, that k has degree-one. Since k is 
analytic, its periodic part has a convergent Fourier series, and so it can be represented as 
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(8) 



with coefficients kj — '"* "" ^'^ 



r,G 



Note that solutions of (7), if they exist, are not unique: given a solution k{9), then k{9 + (f>) is 
also a solution V</> € M. However when co is irrational, continuous conjugacies are otherwise unique, 
apart from this shift in the origin of 9. 

Lemma 1 ([AdlLPOS]). If k G C^CBjM) solves (7) for an irrational rotation number u, then every 
other continuous solution of (7) for the same invariant circle is of the form k{9 + (j)) for some G S. 

Proof. Indeed, if k is another continuous conjugacy, then, since /c(S) = A;(S) = C, there exists 
a (/> such that, e.g., /c(0) = k{(j)), and by (7), this implies k{nuj) = k{(l) + nuj), Vn G Z. Thus 
k{9) = k{9 + (/)) on a dense set in S, and by continuity they agree everywhere. D 



2.1 Automatic Reducibility 

Under the assumption that there exists an invariant circle with the given u, de la Llave and collab- 
orators [dlLGJVOS, FdlLS09] developed an iterative, quasi-Newton scheme to find the conjugacy 
k. The algorithm is initialized with a guess k such that 

fok-koT^ = e, (9) 

and is guaranteed to converge provided that the error, e, is sufficiently small. The iteration proceeds 
by inserting a corrected A: — )• A: + A into (7) and expanding to give 

f{k{e)) + Df{k{e))A{e) - k{e + w) - A{e + w) = o{a^). 

Neglecting the second-order terms and using (9) gives the iterative equation 

A{e + cj)-Df{k{e))A{e) = e{e), (lo) 

that can be viewed as determining A. The resulting function k+A is then an approximate conjugacy 
in the sense of satisfying (9) with a new, presumably smaller, error e. 

Computing A by a direct inversion of the cohomology operator on the left hand side of (10) 
is numerically expensive. However, this linear operator can be partially diagonalized through a 
process called automatic reducibility by [dlLGJVOS]. The idea is that there exists a change of 
variables A(^) = M{9)w{9) where M{6) is a 2 x 2 orthogonal matrix, and such that for the new 
vector w{9), (10) takes the form 

w{e + w) - u{e)w{e) = M^{e + oo)e{e) = e{e), (ii) 

where U{9) is a special upper-triangular matrix. To find U and M, one must solve the matrix 
system 

Df{k{e))M{e) = M{e + uj)u{e). (12) 

To solve (12) the columns of M are chosen to be tangent and normal vector fields of the 
(approximate) circle C, recall Fig. 1. Note that if k solves (7), then differentiation implies 

Df{k{e))Dk{e) = Dk{9 + uj), (13) 

which is the statement that the tangent vector field, Dk, to C is invariant under /. Since the function 
k in (10) will never be an exact conjugacy, Dk will only approximately solve (13). Nevertheless, 
we may use (13) in the Newton iteration (10) incurring only error at second order, see [dlLGJVOS, 
Prop. 7]. The matrix M is now chosen to be 

M{e) = {u{e) vie)) 

where u and v are unit tangent and normal vector fields 

"<^' - P^°'(^'- (14) 

v{d) = Ju{9), 



and J is the standard two-dimensional Poisson matrix. 
With this choice, (12) becomes 



Df{k{e))M{e) = {Df{k{e))u{e) Df{k{e))v{e)) = {u{e + u) v{e + u))(\ ""^^A, (is) 



which must be solved for the ofF-diagonal term of U , the function a{9). Approximate invariance of 
the tangent vector field, (13), implies that the first column of (15) is an identity, at least to second 
order. The second column gives 

a{e)u{e + u) + v{e + w) = Df{k{9))v{e). (i6) 

Taking the inner product of this equation with u{9 + oj) and using orthonormality gives 

a{e) = u{9 + ujfDf{k{e))v{e), (17) 

which determines a. Similarly, taking the inner project of (16) with v{9+u}) and using the definition 
(14) gives the consistency condition 

l = v{9 + ujfDf{k{9))v{9) = {Df{k{9)fju{9 + uj)f v{9) 
= {JDf-\k{9))u{9 + u})f v{9) = {Ju{9)fv{9) = 1, 

where we have used the symplectic condition Df J = JDf^^ and the inverse of (13). Simply put, 
this condition states that the projection of the sheared image Dfv onto v is exactly the normal 
vector V, since area is preserved, see Fig. 1. 

The two rows of (11) now yield skew coupled equations for the components of the vector w. 

wii9 + uj) - wi{9) = eii9) + a{9)w2{9), (18) 

W2(,9 + U) - W2{9) = 62(9), (19) 

where a is defined by (17), and e by (11). 

These two equations can be solved easily in Fourier space. Indeed, each is of the form of a 
cohomology equation, 

w o T^ — w = e, 

that is diagonalized by Fourier transformation. Indeed, when e is analytic and to is Diophantine 
(31), then w is analytic [Mos66] and its Fourier coefficients are 

provided that e satisfies the solvability condition 

eo= /" e{9)d9 = ^, 
Jo 

i.e., that its average vanish. Since the average, wq, is in the kernel of the cohomology operator, it 
can be chosen freely. 



2.2 Summary of the Algorithm 

Beginning with a guess for the conjugacy A;, we compute the error e from (9), and the vector fields 
u and V from (14). Now the transformed error e can be computed from (11) and a from (17). At 
this point the cohomology equation (19) can be solved for W2 using (20), under the assumption that 
the average of 62 vanishes. Even though this assumption is not generally true for an approximate 
conjugacy, ignoring it does not interfere with convergence of the method; indeed as e — )• 0, so does 
the error induced by this inconsistency [FdlLS09]. The arbitrary coefficient W20 can be selected 
so that the average of the right hand side of (18) is zero. At this point (18) is consistent and can 
be solved for wi. The average of wi is arbitrary; indeed, since this gives a contribution u{9)wiq 
to k, it contributes to a shift along the circle, which corresponds to the non-uniqueness of the 
solution, recall Lem. 1. We set iDio = for simplicity. Once w is found, k can be updated using 
k — )• k + M{9)w{6). This process is repeated until the L^-norm of the error (9) is satisfactorily low. 
This algorithm is generally robust and invariant circles can be computed for moderate values 
of £ and many Diophantine u. In some cases the algorithm fails well before £cr, even upon using 
extrapolation in e to update the initial guess. This failure is often the result of aliasing in the Fourier 
spectrum. Aliasing occurs whenever one computes the discrete Fourier transform of a nonlinear 
function, such as g o ki^ of a finite series approximation to k. This error can be ameliorated by 
application of an anti-aliasing filter [TreOO] . We use the quadratic filter 

f kj I •! 7 

'^i ~^ p _ j2 ' l-^l > -^ 

to scale the amplitudes of the Fourier coefficients beyond a threshold J. For the twist case, rii(y), 
we set J = 2N for N Fourier modes, but we set J = ^N for ^2 and other nonlinear frequency maps 
to help compensate for the additional nonlinearity. This filter is applied to k once per iteration, at 
the beginning of the Newton step, prior to the computation of the error e with (9). 

To compute invariant circles for large e, we use continuation from e = 0, where k is trivial. We 
begin with the increment Ae = 0.01 and with N = 2^ Fourier modes. The algorithm is deemed 
to converge when the L^-norm of the error e is less than the tolerance 10~^^. An initial guess for 
the conjugacy at e -|- Ae is obtained by linear extrapolation. This is continued until the algorithm 
fails to converge within the specified tolerance. Upon this first failure, the step size is reduced to 
Ae = 0.0005 and the number of Fourier modes is doubled. At each successive failure the number 
of Fourier modes is again doubled; however, we found it significantly faster and more accurate to 
keep the step size constant. For reasonable accuracy and speed, we let the algorithm exit when 
the number of Fourier modes exceeds 2^^ (conjugacies with millions of Fourier modes have been 
computed [AdlLPOS, OP08]). This failure occurs below the critical value, £cr, at which the invariant 
circle is first destroyed. Indeed, a comparison with Greene's criterion (see App. B) for Chirikov's 
map, as we discuss in §3, shows that convergence typically fails when e r--^ £ct — W~^- 

3 Detecting Critical Circles 

The quasi-Newton scheme discussed above provides a method to compute an analytic invariant 
circle when it exists. We will use the scheme of §2 to estimate the value £cri^) such that there is 



an invariant circle with rotation number u for all < e < £cri^)- As in Conj. 1, we expect that 
when £ exceeds £cr the invariant circle is replaced by a cantorus. 

Since the conjugacy loses continuity at £cr{^), this transition can be detected by the blowup of 
a Sobolev norm [AdlLP05, OP08, CdlLlOa, CdlLlOb]. To detect the transition, we use a related 
seminorm defined as the L^-norm of the m derivative: 

\\k\t^\\D^k\\l = Y,{2n\j\f^\k\^ (21) 

j 

where k are the Fourier coefficients of A;. As e — )• £cr, \\k\\m — ^ oo, and we assume it does so 
asymptotically as 



(22) 



To compute the three parameters in the asymptotic form (22), we use three consecutive (e, ||A;||m) 
pairs from the Ae steps of the continuation method. One reason for using a fixed Ae, as described 
in §2.2, is that we find the estimate of the pole position to be more accurate than if a variable 
step size were used. We typically compute (21) for the angle component of k, though the action 
component gave similar results. 

Numerical computations for the rotation number 7"^ and several examples of (1) indicate that 
(22) is a good fit for m = 2 when £cr — e ^ 10~^, see Fig. 2(a). More generally we find that (22) 
applies for other values of m with b Ri m — 1. 

The choice of m significantly affects the error in the approximation of £cr- In Fig. 2(b) we 
compare the estimates of £cr from (22) for Chirikov's map with those obtained numerically from 
Greene's criterion (recall App. B). For the invariant circles with rotation numbers 7 and \/2, and 
using orbits of periods up to 30, 000, Greene's method gives 

£,,(7) = 0.9716353(1), ^ ^ 

(23 
£^riV2) = 0.957447(6), 

where the error (the number in parentheses is the uncertainty in the last digit) is estimated from 
the extrapolation.'^ We observe that the smallest numerical error for the seminorm fit occurs with 
the choice m = 2. Fixing the maximal number of Fourier modes to 2^'^, mainly for computational 
speed, we typically obtain e^t with an error less than 10~^ (more details will be given in §5). 

When the invariant circle is not symmetric, we cannot compare with the results from Greene's 
residue criterion (a periodic orbit finder would require a two-dimensional search and be much less 
well-behaved). Nevertheless, we select m = 2 for the computations in the rest of this paper. 
We estimate the error in the calculation of £cr as the difference between the last two, three-step 
approximations obtained with the step size Ae = 0.0005. Note, however, that this occasionally 
underestimates the true error in £cr by a factor of 2 — 8. 

4 Conjugacies of Near- Critical Circles 

In this section we compute the conjugacy (7) for circles with golden mean rotation number for a 
variety of maps of the form (1). We use the seminorm technique to determine when these circles 



'Quadruple precision computations give Ecri'y) = 0.9716354063(2) [Mac93, §4.4] 
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Figure 2: (a) Blowup of the seminorm ||fca;||2 for the uj — 7~^ invariant circle for three cases of (1). The upper 
(blue) curve corresponds to Chirikov's standard map. The lower two curves have the force (25) with %j) — 7r/4. The 
middle (red) curve is for the frequency map (3) and the bottom (green) is for (5) with 5 = 0.3. In each case the 
horizontal axis is logarithmic based on the best estimate of e^r from (22). (b) Error in the pole location for fits 
to (22) for seminorms with varying m for Chirikov's map for two rotation numbers, compared to Scr from Greene's 
criterion. For both rotation numbers, the smallest error occurs with m — 2. 

are destroyed, comparing to Greene's criterion when the map is reversible, but also generalizing to 
maps that do not have twist or reversing symmetries. 

We begin by computing the embedding of the oft-studied golden mean invariant circle for 
Chirikov's map, i.e., for g = gi, (2), and fi = Qi, (3). The algorithm of §2 converges up to 
£ = 0.9695, which from (23) is £cril) — 0.0021. The resulting embedding at this parameter value, 
as shown in Fig. 3, has error e{6), (9), with L2-norm of (10)~^^ and L°°-norm of 1.7(10)~^^. 

Since Chirikov's map has twist, Aubry and Mather theory implies that each invariant circle 
becomes a cantorus upon destruction. Moreover, these sets are ordered in the sense that 9 i— )• kx{0) 
is a monotone circle map. The implication is that upon destruction of the circle, the extension of 
the conjugacy to the Cantor set is the inverse of a devil's staircase; consequently, kx{9) — 9 must 
develop positive jump discontinuities at a dense set of 9 values corresponding to the location of 
the gaps in the cantorus. Although the Fourier method does not converge for the super-critical 
case, the beginning of this metamorphosis can be seen as e — )• £cr- Indeed the largest, incipient 
discontinuity seen in Fig. 3(a) occurs at 9 = ^, which, since the computations give kx{^) ~ ^, 



also corresponds to x 



The remaining localized regions with large slope occur along the orbit 



?t = I + tUJ. 



The locations of the gaps can be explained by recalling that an invariant circle corresponds 
to a minimizing state of the Frenkel-Kontorova energy (4). The potential for this map, V{x) = 
— ^;^cos{2iTx), has a single maximum at a; = j! thus, the AI theory implies that the largest gap 
should form around this maximum, and its images form a bi-infinite family (one "hole") that 
corresponds to the gaps in the cantorus. The first few forward and backward iterates in this family 
are indicated by the vertical lines in Fig. 3(a). Continuing this process to larger iterates appears 
to account for all the local peaks in the derivative of kx; thus it appears that this cantorus has one 
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Figure 3: Components of the conjugacy for the oj — j~^ invariant circle with of Chirikov's standard map from a 
computation with 2^^ Fourier modes, (a) The function kx{6) — 9 and the first three forward and backward images 
oi 9 = |, indicated by the sohd (red) and dotted (green) Unes, respectively, (b) The action component ky(9) is even 
about uj/2. 

hole [MMP84], a fact that can be proven near the AI hmit [BM94]. 

The standard map has two independent famihes of reversors, see App. C: it is doubly-reversible. 
The first reversor arises from the oddness of gi and the second from the oddness of ili. As is shown 
in Cor. 4 in the appendix, the first reversor imphes that the function kx{0) is odd about some point 
ip; using (42), we computed (p = —kxo ~ —8.130(10)^'^. Since tp is so smah, the graph in Fig. 3(a) 
looks like it is odd about 0. Corollary 4 also implies that the action ky{6) is even about ip + gW, as 
is visually apparent in Fig. 3b. More precisely, we found that the identities (42) are satisfied up to 
an error with L°°-norm of 9.9(10)"^^, comparable to the accuracy of the computation of k itself. 

A second, doubly-reversible version of (1) is the generalized, two-harmonic twist map with the 
force 

g2{x; "0) = 5^ (sin(^) sin(27rx) + cos(0) sin(47rx)) , (24) 

keeping the frequency map, $7i. The breakup of the golden circle for this map was studied first by 
Greene et al [GJSS87], and later the formation of cantori and turnstiles were studied in [KM89, 
KM94, BM93, BM94, LC06].^ This two-harmonic map reduces to Chhikov's map when V = f • The 
two maps are also equivalent when ip = or ir under the rescaling transformation (x,y) i— )• 2{x,y), 
so that the parameter and rotation number are effectively doubled. Finally, since g{x + ^, —tp) = 
g{x.,il))., it is sufficient to study the parameter range ^ G [0,'/r]. 

Two examples of near-critical invariant circles for the two-harmonic twist map are shown in 
Fig. 4. Since this map has the same reversors as Chirikov's map, kx is still odd (since kx{^) ~ 0, 
essentially about = 0). When f/; is small, as in Fig. 4(a), it appears that there is a symmetric 
pair of two equally large incipient gaps (near 9i = 0.3090 and 0.6910), while for larger -0, Fig. 4(b), 
there is a single, largest gap near 6 = 0.5. In both cases, it appears that there is a single orbit, 

^In these papers, the parameterization fci = esin(?/)) and fe = —2ecos{ip) was used for eg2- 
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g,(x,0.75ji),Q,(y)J 
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Figure 4: Plots of kx{6) — 6 for near-critical invariant circles for (1) with (.3) and (24) and uj = j~^ ■ Also shown are 
the first three forward (red, solid) and backward (green, dashed) images of the largest gap. Insets show the shape of 
the potential V{x) for the two values of tp. 



Oi + toj, of gaps: for example, two largest peaks in Fig. 4(a) are separated by one iterate. Thus the 
resulting cantori will have one "hole." 

The difference between the two cases is correlated to the structure of the potential V{x) = 
J 92dx — see the insets in Fig. 4 — which has critical points at x = 0, ^ and any solutions of tan-;/' = 
— 2cos(27r2;). When < V' < arctan2 ss O.SSvr, V has two wells, the deepest at x = 0, separated by 
maxima of equal height; this is the case for Fig. 4(a). The maxima collide at tan'0 = 2 annihilating 
the well at X = 2- For | tau'^l >2,V has single well. Finally, when 0.657r ~ vr — arctan2 < ip < tt 
the potential has two equally deep wells separated by differing maxima, the largest at x = ^j as 
shown in Fig. 4(b). The AI limit, e — t- oo, for a double well potential has infinitely many cantori, 
selected by the fraction of points that lie in each well [KM94, BM94]. For these cases, there are 
cantori with two families of gaps; however, the second hole is formed beyond the initial breakup of 
the invariant circle. Note that when tan'0 = ±2, the potential has degenerate critical points and 
persistence of AI states with these points occupied cannot be guaranteed. 

When the maxima of the potential collide at tan^ = 2, we expected that the symmetric pair 
of largest gaps shown in Fig. 4(a) would also merge. However Fig. 5 shows that this merger does 
not occur until ip ~ 0.46-7r, well after the collision. Note that as the largest gaps coalesce, they still 
appear to lie on a single orbit, though the number of iterates between them changes with ip. For 
example, for ip = 0.447r, in Fig. 5(a), the largest gaps are two iterates apart, but for ip = 0.467r, in 
Fig. 5(c), they are five iterates apart. 

The set of critical parameters {Scrjip) for the golden circle of this map, computed using the 
seminorm method, is shown in Fig. 6(a). For -(/; = f , ^cr ~ 0.971635 since g2ix, f ) = gi{x). In 
addition Scr ~ 0.40236 when ^ = or vr — this is half the critical parameter for the uj = ^ circle 
of Chirikov's map [GJSS87]. When ip G [0, |], the critical set exhibits a Cantor set of cusps as 
was first observed using Greene's criterion in [KM89]. The cusps are related to a complex set of 
symmetry breaking bifurcations of periodic orbits that limit onto the invariant circle [KM94] . 
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Figure 5: Near-critical conjugacies, kx{9) — 6, for the 7~^ circle of (1) with (72 and f7i, and the values of ip and 
e indicated. Insets show the potential V{x) for these parameter values. Vertical lines show forward (solid) and 
backward (dashed) iterates of the largest gap. 

To begin an exploration of the breakup of circles in maps with less symmetry, consider a two- 
harmonic force similar to (24), but with a phase shift so that it is no longer odd: 



5(3 (x; ip) = 2^(sin(-(/;) sin(27rj;) + cos{ijj) cos(47rx)). 



(25) 



With the odd frequency map Oi, the map still has one reversor S2, (38), which conjugates the two 
circles with rotation numbers ibu;. As before, 53(x + ^j—ip) = gs{x,ip), so the parameter ip can 
be restricted to the range [0, vr]. Moreover, since g^{—x,'K — -0) = —gs{x,'ip), the transformation 
(x, y, ip) — )• {—X, —y, IT — ip) leaves the dynamics invariant, so we need consider only ip £ [0, |]. 

Two examples of near-critical invariant circles for this map are shown in Fig. 7. Since the force 
(25) is not odd, the map (1) no longer has the reversor (37), and the resulting invariant circles 
will not be invariant under Si. Indeed, the conjugacies shown in the figure no longer show this 
symmetry and the relations (42) no longer hold. 
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Figure 6: Critical parameter set for the golden circle of (1) for the frequency maps and two-harmonic forces shown, 
computed using the seminorm. The axes are the amplitudes of the two Fourier modes of eg. The dashed lines 
correspond to degenerate AI limits and for (a) to ^ = 0.467r where the two symmetric gaps coalesce in Fig. 5. Insets 
show representative V(x). 





ft(x,0.25;t), n,0') 
8=0.318 




Figure 7: Angle component of the conjugacy for a near-critical circle of the generalized standard map with gz and 
f^i, for uj — 7~^ and two values of i/i. These embeddings do not exhibit the symmetry seen for the reversible case, 
but do display an orbit of incipient gaps. 

The potential for gj, has two wells when ?/; € [0, ^) and a single well for V' £ [f , tt]. In both panes 
of the figure, the potential has a single well and the largest gap in kx is formed near the maximum 
of V . For ^ < ^, the next-largest gap is formed near the smaller maximum of V . This gap can 
still be seen near = 0.2 for the degenerate case of Fig. 7(a), and the corresponding x = kx{6) is 
near the degenerate critical point of y at x = 0.25. This gap is the foward image of the largest 
gap and is much larger the remaining gaps. This secondary gap is still present in Fig. 7(b), but is 
significantly smaller, giving some evidence for the influence of the smaller maximum on the size of 
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the gaps. 

The critical set for the golden circle of this map is shown in Fig. 6(b) (the black curve) for the 
range ip £ [0, |]; this set is symmetric under reflections about the horizontal and vertical parameter 
axes. Just as for Fig. 6(a), the points at "0 = and V' = f correspond to Chirikov's map. However, 
this critical curve no longer has a Cantor set of cusps — indeed, it seems to be smooth apart from 
cusps at ^ = and n, corresponding to the doubled Chirikov map. 

We now consider several maps with the frequency map Q2, (5), for which twist reverses at y = 0. 
For this case, there are typically two circles for each rotation number in the range of i^2{y, S)- When 
u) ^ —5 one circle is contained in the positive twist region, y > 0, and one in the negative twist 
region, y < 0. The most interesting circles, however, cross y = so that the twist condition is 
locally violated. The breakup of these shearless circles for the standard nontwist map (with gi) 
has been much studied [AWM03, AdlLP05, FWAM06, WFPll], so we do not examine this map in 
detail here. 

Instead we consider a nontwist map with force (73 . It would seem that destroying the oddness of 
both g and O, would eliminate both reversors of (1). However, since gs is even about j, g3{^+x, tp) = 
g3{^—x,ip), and ^2 is even about y = 0, this map has a reversor, the map 5*3 of (39). This reversor 
conjugates the positive and negative twist circles; in particular they have the same Ecr- Moreover, 
the transformation (x, y, ^) — )• (x, — y, tt — ip) leaves the dynamics invariant, and combining this 
with the reversor allows us to restrict ip £ [0, |]. 

Invariant circles of this map with oj ^ —6 are very similar to those of a twist map, see for 
example the positive twist circle in Fig. 8. The conjugacy is very close to that for the twist case, 
O = Qi, shown in Fig. 7. Indeed, since {ky{6)) w 0.9581 and varies by only ±0.023, the twist 
is nearly linear over the range of oscillation of y. Moreover, since DVt2{y,5) = 2y ~ 1.9161, the 
twist is almost a factor of two larger than that for f^i — this accounts for the fact that Ect has been 
decreased by factor very close to two from the ili case. 








g,{xfi.25Ti), n2Cv,0.3) 
8=0.164 



Figure 8: Angle component of the conjugacy for near-critical circle of the generalized standard map with g^ and 
^2, for uj = 7~^ and y > 0. 

Circles that experience a stronger variation in twist display significantly different behavior than 
those of the twist case. This is especially true of those circles that cross the y = axis, where 
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the twist vanishes. Figure 9 compares circles for Qi and ^2 with the same noble rotation number 
w ~ —0.2793. Both still have one orbit of gaps; however for Q2, the dominant gap is no longer 
associated with the largest maximum of V. Thus, though it appears that the invariant circle will 
become a "one hole" cantorus upon destruction, traditional anti-integrable theory may not provide 
insight into its formation. Note also that Scr for the nontwist case is 4.02 times larger than that for 
the twist map, no doubt related to the fact that (Dri2(y, 0.3)) = 2kyQ = 0.2083 is 4.8 times smaller 
than DQi = 1. 



I 




g3(x,0.25;t), n,(y) 
8=0.351 




g,{x,0.25n), 0,0,0.3) 
8=1.412 



Figure 9: The conjugacy k^ for a near-critical invariant circle for g3{x, j) with w — —(12 + 197)/(43 + 687) and 
r^i (a) and f22 (b). The invariant circle for the nontwist case (b), crosses the y-axis. 

Finally, we explored a map that we conjecture has no reversors. For this map we use the 
frequency map introduced by Rannou [Ran74] 

^3{y) = -y+ ^(cos(27ry) - 1), (26) 

that is neither even nor odd. This frequency map, like Q2, does not satisfy the twist condition. 
Rannou used this in conjunction with an even force that that has a nonzero integral so that the 
map has nonzero net flux and no rotational invariant circles when e 7^ 0. To permit invariant 
circles, but to eliminate additional symmetry, we use instead a force that is neither even nor odd, 
has a zero average, and a full Fourier spectrum: 

54(x; tp) = -z— sin (27rx + ip cos(47rx)) . (27) 

The generalized standard map with (^3 and 54 is, to the best of our knowledge, completely nonre- 
versible and has only a trivial symmetry group. 

An example of a near-critical golden circle for this map is shown in Fig. 10. This circle also 
appears to have one hole; the principal gap forms near 6 = 0.44, close to the (single) maximum of 
the potential, as shown in the inset. 

In conclusion, for all of the cases that we investigated, the near-critical invariant circles show 
the formation of a single family of gaps, the largest of which is usually associated with the global 
maximum of V{x). This presumably leads, as in Conj. 1, to the formation of a cantorus for e > Scr, 
even when the twist condition is not satisfied. 
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Figure 10: Angle (a) and action (b) components of the conjugacy for a near-critical invariant circle with uj — j ^ 
of a generalized standard map with no known symmetries. 

5 Robustness of Noble Invariant Circles 

An invariant circle is locally most robust in a one-parameter family /^ if it exists for larger e 
than any circle in some neighborhood in phase space. The Greene-MacKay renormalization theory 
[Mac93] leads to the conjecture that the locally most robust circles for twist maps have "noble" 
rotation numbers, that is 

a;(M,7) = ^, u = {1,^)M, (28) 



for a unimodular matrix M. In this case z/ is an integral basis for the ring ^(7), and the collection 
of nobles is projectively equivalent to the set of integral bases of ^(7), recall App. A. 

It is not unexpected that the nobles are robust; for example, the size of the perturbation 
allowed in KAM theory is proportional to the Diophantine constant c(w) (31), and any noble has 
the maximal value 0(7) = 5~^'^. The most conclusive evidence for the robustness of the nobles is 
the numerical study of MacKay and Stark [MS92] for the standard map. 

Here we investigate the relative robustness of nobles relative to other quadratic fields for the 
various generalized standard maps introduced in §4. The six real rings Z(r) with the smallest 
discriminants D{r), or equivalently the largest Diophantine constants, are shown in Table 1. To 
characterize the computational accuracy, we first compute the critical parameter value of Chirikov's 
map for circles in each ring using both the residue and the seminorm methods. Table 1 shows that 
the golden mean is, of course, the clear winner with the largest critical value, (23). 

The error in each computation of Scr was estimated by comparing the extrapolated values from 
the last two iteration steps. The computations indicate that Ecr is about an order of magnitude 
less accurate using the seminorm method; however, the error is occasionally overestimated by this 
calculation as discussed in §3. It is not surprising that error for both methods tends to increase 
with discriminant: circles from fields with larger D tend to be closer to low-order resonances. For 
the residue method, each periodic orbit is nearer to those with lower period, making the orbits 
more difficult to compute, and for the conjugacy method, a larger discriminant implies that the 
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[3,1] 



e^Jr) 



0.9716358(1) 

0.957447(6) 

0.87608(1) 

0.89086(2) 

0.91573762(4) 

0.77242(5) 



.(r) 



0.971638(1) 

0.95744(7) 

0.8756(6) 

0.8905(3) 

0.9143(3) 

0.7720(8) 



{As 



cr/rms 



1.6(10)-'^ 
3.0(10)-^ 
4.3(10)-^ 
3.3(10)-3 

1.0(10)-^ 
2.1(10)-^ 



Table 1: The six quadratic fields Q(r) with the smallest discriminants, D. Each is generated by the irrational r 
with the periodic continued fraction (CF) shown. The 4*** and S**" columns give ecr{r) for Chirikov's map computed 
by the residue method (for orbits up to period 30, 000) and the seminorm method, respectively. The parentheses 
indicate the estimated extrapolation error in last digit. The last column is the root mean square difference between 
the two computations for 256 rotation numbers in each field. 



denominators in the solution of the cohomology equation (20) will be smaller, inducing larger errors 
in the Fourier coefficients. 

To study local robustness, we construct a set of nearby rotation numbers from each ring, using 
the 256 rationals in [0, g] from the Farey tree up to level eight, see App. A. Each pair of neighboring 
rationals, {pi/qi,P2/Q2), gives a basis vector u = (pi +P2r, qi + q2i") for each of the six rings. The six 
corresponding frequency ratios u}{M,r) (28) lie the interval between the neighbors. These interval 
widths range from 0.0002 to 0.1 and define the "local" neighborhoods for the robustness test. The 
last column of Table 1 shows the rms difference between the residue and seminorm computations 
for the 256 irrationals in each ring for Chirikov's map. 

For a given pair of neighboring rationals, we compare the values of Ecr for each of the six rings 
to determine which is locally most robust. Since the computations of Ecr sue uncertain, we are 
unable to determine a winner when the difference between the most robust and other circles is too 
small. The approximation to Ecr combined with the estimated error generates a critical interval, in 
which the true Ecr likely resides. A circle is deemed to be locally most robust if its entire critical 
interval is above the critical intervals of the other five circles. If the most robust interval overlaps 
with any other interval, we cannot conclusively declare a winner. We also discounted any result for 
which the width of the critical interval for the most robust circle was greater than 10~^. 

Table 2 shows the number of times each ring was most robust in a comparison with the other 
rings. For the residue computation the nobles are most robust 69% of the time, and for the 
seminorm method, 67%. Errors caused 2-3% of the trials to be discarded. These results are similar 
to those of MacKay and Stark [MS92]. There are cases in which numbers in Q(-v/2) appear more 
robust than a noble, but MacKay and Stark found that when the noble at a particular level on 
the Farey tree was not most robust, there was always a noble at a deeper level that was. Thus the 
cases in which a noble does not win the contest appear to result from the neighborhood being too 
large. 

If we remove the field Q(7), then a comparison of the remaining five fields (the second and 
seventh rows of the table) shows that Q(v2) is most robust. Continuing this comparison for the 
remaining fields shows that robustness seems to be a monotone function of the field discriminant. 
However, for the seminorm method, the uncertainties are too large to determine a winner amongst 
the four largest discriminant fields. 
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Table 2: The relative robustness of invariant circles of Chirikov's standard map with rotation numbers in different 
algebraic rings. The numbers in each column indicate the number of times a circle with rotation number in that ring 
was most robust for 256 trials. The final column is number of trials for which a winner could not be conclusively 
determined. 

It has long been conjectured that noble invariant circles should be locally most robust for maps 
with more general forces, with nonmontone frequency maps, and without reversing symmetries. 
We investigated this using the seminorm method for the maps of §4 and found that in each case the 
nobles appeared to be locally most robust to nearly the same degree that they are for Chirikov's 
map. For example. Table 3 shows the results for the generalized standard map with g^ and ili for 
i/j G [0, 7r/2]. Regardless of the choice of ip, the noble circles were once again locally most robust 
at least 60% of the time. We suspect that — ^just like for Chirikov's map — whenever one of the six 
non-nobles was a winner, there is a noble below it on the Farey tree that wins, but we did not 
examine this in detail. 

As a second example. Table 4 shows the robustness data for the nontwist map with frequency 
map 0.2- In this case the symmetries of the map that allow us to restrict to a; G (0, ^) are gone. 
We therefore studied the 512 rotation numbers generated from the rationals on level 9 of the Farey 
tree in the interval (0, 1]. We once again found that circles with noble rotation numbers are locally 
most robust at least 55% of the time. For this case the errors were larger, and up to 9% of the 
comparisons had to be discarded. 

In conclusion, our results support Conj. 2, that invariant circles with noble rotation numbers 
are generically, locally most robust for the generalized standard map (1). 

6 The Critical Function 

The destruction of the last rotational invariant circle is a significant event in the dynamics of area- 
preserving maps since this corresponds to the transition to global chaos and unbounded orbits. 
In this section we will identify which of the nobles is globally most robust for the some of the 
generalized standard maps studied in §4. Since, as we confirmed in §5, the noble circles appear to 
be locally most robust, we will confine our search to these rotation numbers. 

The last invariant circle is the global maximum of the critical function Ecri^^)- The conjecture 
that the golden circle is most robust for Chirikov's map has been supported by many computations 
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Table 3: Relative robustness of circles from six quadratic fields for the generalized standard map with g3{x,ip), 
Oi(3/) and 10 values of il) using the seminorm method. The number in each column corresponds to the number of 
times each of the six rings of Table 1 was most robust from 256 trials. The final column is number of trials for which 
a winner could not be conclusively determined. 
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Table 4: Relative robustness of circles from six quadratic fields for (1) with gi{x,ij)) and 572(1/, 0.3) and 10 values 
of ?/; using the seminorm method. The number in each column corresponds to the number of times each of the six 
rings of Table 1 was most robust from 512 trials. The final column is number of trials for which a winner could not 
be conclusively determined. 
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using periodic orbits [SB82, BPV90, Mur91]. It is also supported using the seminorm technique, as 
can be seen in Fig. 11(a). Here we computed ecr{^) for 256 noble frequencies in (0, g) using level 8 
of the Farey tree as in §5. The resulting picture is indistinguishable from that obtained using the 
residue method and has maximum value ecr(2 — 7) = Scri'j) of (23).^ 
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Figure 11: Critical functions for (a) Chirikov's standard map and (b) the standard nontwist map for S = 0.3, using 
the seminorm method. 

Determining the globally most robust circle for the standard nontwist map, with O2 and gi, is 
complicated because of the lack of periodicity in frequency space. Figure 11(b) shows the critical 
function for —6 < to < 1 for 512 nobles generated from level 9 of the Farey tree with root [j, j] and 
112 negative nobles greater than —6. When w > 0, this critical function has a similar structure 
to Chirikov's map; however — as has been often observed [MorOO, SCL^09] — the circles become 
increasingly more robust as w — >• —6. Indeed we found the most robust circle to be that with the 
smallest co. Unfortunately, the extrapolation errors for the estimate of £cr tend to become larger 
in this limit and we discarded the estimates for 44 circles for which the error was larger than 10~^. 
A similar picture of the critical function was obtained for the nontwist map with the force (73: the 
most robust circle was always found to be the one with rotation number closest to —6. 

The critical function is considerably more complex for the multiharmonic twist map, with Qi 
and §2- Figure 12(a) shows that the graph of £cr{^) resembles that for Chirikov's map when 
ip G [f )7r]; however, when tp G [0, |], Ecri^) is much less regular, see Fig. 12(b). Recall that this 
region of parameters corresponds to the Cantor set of cusps in the critical set of the golden circle, 
e.g.. Fig. 6(a). We also found that £crii^) is especially erratic when two gaps collide, as we observed 



near ip 



0.47 in Fig. 5. 



^Symmetries of Chirikov's map, recall App. C, imply that £cr{(^) = £cr{n ± i^) for n G Z, so it is sufficient to 
consider uj £ [0, |]. 
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Figure 12: Critical function for the generalized standard map with Qi and (a) 52(2^, 2.3484), or (b) 32(2;, 0.7776) 
for 256 noble rotation numbers. 

7 Most Robust Circle 

To provide evidence for Conj. 3, we studied the rotation number of the most robust circle 



Wn 



argmax{ecr(^)} 



for two parameter families of generalized standard maps. In particular, we studied the dependence 
of ujmax on the parameter ^ in our various force models. To do this, we first obtained a rough 
estimate for the interval of rotation numbers that contains the most robust circles over a range of 
■0. Then, for each value of ip, we select the most robust circle from the 256 nobles on level 8 of the 
Farey tree whose root corresponds to this interval. 

For example, for the frequency map Qi, periodicity allows us to restrict u to the interval [0, j]- 
The dependence of comax on the parameter ip for the force g^ is is shown in Fig. 13(a). Surprisingly, 
of the 256 nobles examined, only five are most robust over most of the range of tjj; their rotation 
numbers are shown in the figure. The critical e for each of these circles varies smoothly with ^, as 
shown in Fig. 13(b). 

The values of Umax {'4') shown in Fig. 13(a) follows a predictable pattern. As noted in §4, this 
map reduces to Chirikov's map at "0 = |, and to the same map on half the spatial scale at ^ = 0. 
Thus the circle with rotation number 2 — 7, whose continued fraction expansion is [0, 2, 1], is most 
robust at 77. When ■0 = the critical function has equal global maxima at the "half-noble" rotation 



2- 



numbers ^('^ i 7)- For the range [0, j] there are maxima at 



i(7 
i(2- 



-1) 

7) 



[0,3,4], 
[0,5,4]. 



(29) 



For moderate 0, the observed values oicomax move logically from 2—7 towards one of the half-nobles, 
in the sense that their continued fraction expansions become increasingly similar to those of (29). 
This behavior persists for small "0, as shown in the vicinity of 5(7 ~ 1) ^^ Fig. 14. The evolution is 
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Figure 13: Globally most robust circles for ili and g3{x,tp) for 100 values of t/j G [0, ^]. Labels show the continued 
fraction expansions of the rotation numbers for five most robust circles, (a) The rotation number of the most robust 
circle, (b) The critical sets, {(cosi/),sin7/))£cr(^)} for the five most robust circles of (a). 

not inonotonic — indeed, the rotation number of the most robust circle appears to oscillate about 
2 (7 — 1). For tp < 10"'^ the rotation number of the most robust circle in the figure is [0, 3, 4, 4, 2, T], 
the noble closest to 9(7 ~ 1) amongst the 256 in the sample. 



[0.3,4,4,2,1] 



[0.3.42,1] [o_3_4j, 
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[0.3,4,6,1] 



[0,3,4,3,2,1] 



[0,3,2,2,1] 




Figure 14: Most robust circles for Qi and g3{x,ip) for 50 values oi ip € [10~®, 1] sampled from 256 circles with 
noble rotation numbers in the interval [f , |]. (a) The rotation number of the most robust circle. Labels show the 
continued fraction expansions of the rotation numbers of the most robust circles, (b) The critical curve, ecr(w) for 
the most robust circles of (a). 

Thus we see that the rotation number of the most robust circle seems to remain constant on 
intervals of the parameter ip, at least for the twist map with force 53. We examined several other 
multiharmonic twist maps in order to determine if this behavior was typical. Figure 15 shows 
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^maxiip) for the generalized standard maps with frequency map fJi and forces 54 and 

sin(27r2;) 



95{x;ip) 



2tt{1 + V'sin(47ra;))' 



(30) 



Since both maps reduce to Chirikov's map when -0 = 0, the golden circle is most robust at this 
point. As ip grows, the rotation number of the most robust circle grows monotonically in steps 
approaching 2- It is interesting that the rotation numbers of the most robust circles for the two 
maps are identical, though the range of ip is shifted. It seems surprising that the circles with 
rotation numbers near ^ are most robust as V' ~^ 1- A partial explanation for this is that the 
Fourier expansions of both g^ and 35 have only odd modes, so in particular the (1, 2) resonance is 
suppressed. 




Figure 15: (a) ujmaxXi') and (b) Ecriij^max) amongst 256 nobles in [0, |] for the twist map with forcing g4,{x, ip) (27) 
(solid black) for xp G [0, 1] and 35(0;, tp) (30) (red open circles and dashed line) for 7/; G [0, 0.7]. 

The behavior of the most robust circles for the fully nonreversible case, with il = (^3 and g = §4, 
is similar to the twist case, see Fig. 16. For this map, the globally most robust circle always had a 
small rotation number and the proximity of these circles to resonance led to larger errors. To focus 
on the robust region, we restricted to the 64 nobles from level 6 of the Farey tree with the root 
[0, ^], and to improve the accuracy, we used up to 2^^ Fourier modes and smaller step sizes in e, 
beginning with Ae = .005, decreasing to As = .0001 when that step size failed. The results indicate 
that ujmax is again constant over intervals of ^p. Also, just as for the previous cases, Scri^max) is a 
monotone decreasing function of ip. 

In each of the maps that we examined so far, ecr('^jV') is a smooth function of ip for each 
fixed noble rotation number. The resulting function cvmaxiip) appears to be piecewise constant, 
supporting Conj. 3. However, for the doubly-reversible, two-harmonic map (f^i and (72); £cr(7)'0) 
was seen to be a highly irregular function when ip G [0, ^] (recall Fig. 6). A similar irregularity 
in Ecr also occurs for other rotation numbers, as suggested by the graph of the critical function 
Fig. 12(b). This irregularity also generates erratic behavior in comax for the same range of ip, see 
Fig. 17. Interestingly, when ip £ [f ,7r], uJmax is piecewise constant and ecrit^max) is smooth, as for 
the previous examples. Consequently, the irregular behavior of the most robust circles appears to 
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Figure 16: (a) LJmax(t^>) and (b) Ecrio^max) for the generalized standard map with with Q, = 0,^, g = 34, estimated 
using 256 nobles in [0, |] for each of 100 values of tp £ [0, 1]. 

reflect the same complex set of symmetry breaking bifurcations that are responsible for the Cantor 
set of cusps in the critical set for the golden circle. 

Thus we suggest that the conjectured piecewise constant behavior oiojmax will occur generically 
in families maps that do not have special symmetries that give rise to irregular dependence of Ecr 
on a parameter. 
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Figure 17: (a) uJmax{ip) and (b) Scr{oJmax) for the generalized standard map with f2i and 32 as a function of 1/), 
estimated using 256 nobles in [0, |] for each of 200 values oi ip € [0, tt]. 
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8 Conclusion 

The Fourier method, outhned in §3, provides an accurate and efficient approach to find the embed- 
ding, k(9), for invariant circles of area-preserving maps with a given rotation number. Using the 
blow-up of a Sobolev seminorm, we obtained reasonably accurate estimates for the parameters at 
which a Diophantine invariant circle is first destroyed; however, the errors in this computation are 
larger than that for Greene's residue criterion. Nevertheless, the Fourier method is easily general- 
izable to maps without the reversing symmetry that makes the computation of periodic orbits for 
Greene's method relatively easy. 

In this paper we have provided some support for the conjectures discussed in §1. For example, 
in all of the cases we examined, the near-critical conjugacies show the incipient formation of gaps 
that signal the mechanism for destruction of the circle and the creation of a remnant cantorus. 
This provides some evidence that — even when Aubry-Mather theory does not apply (i.e., the twist 
condition is violated) — cantori form in area-preserving mappings. It would be nice to generalize 
Aubry's anti-integrable limit to such maps to confirm that they have cantori near e = oo. 

In §5 we gave numerical evidence that the relative robustness of a circle depends monotonically 
on the discriminant of the algebraic ring of its rotation number. Thus the noble invariant circles 
appear to be locally most robust as conjectured in Conj. 2. Again, this conjecture has much support 
for reversible twist maps, but the seminorm technique has allowed us to extend it to nonreversible 
and nontwist maps. 

Our investigation was limited to rotation numbers in quadratic algebraic rings. We hypothesize 
that this would extend to more general irrationals, however, to our knowledge, this has never been 
studied. The extension of this relationship to tori in higher dimensional symplectic and volume- 
preserving maps is also an open question. 

Finally, support for Conj. 3 was given in §7. There we observed that — except for the doubly- 
reversible, two-harmonic map — the rotation number of the globally most robust circle seems to be 
piecewise constant when a second parameter of the map is varied. The exception to this seems to 
be related to extra symmetry which gives rise to a highly irregular critical function. 
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Appendices 

A Quadratic Irrationals and the Farey Tree 

The rotation number a; of a persistent invariant circle in the standard version of KAM theory 
satisfies a Diophantine condition, i.e., there is a c > and an s > 1 such that 

q^\quj — p\ > c, (31) 

for all p S Z and (7 G N. Numbers that are Diophantine for s = 1 are called badly approximable, 
and the largest asymptotic value of c for s = 1 is the Diophantine constant: 

c{uj) = lira, inf q\qu} — p\. (32) 

A number cj € C is a quadratic algebraic integer if it is the root of a monic quadratic equation 
with integer coefficients [HW79]. The integer to generates both an algebraic field, Q{uj) = {a + bcv : 
a,b £ Q}, and a corresponding ring of integers, Z(a;) = {a + 6a; : a, 5 G Z}. A basis for this ring 
is referred to as an integral basis. Every real, quadratic integer is Diophantine, and its constant c 
can be easily determined. 

Theorem 2 ([Cas57]). If u: is a real, irrational root of the monic polynomial x^ + bx + c, with 
a, 6 G Z, then 

c{^) = ^, (33) 

where D = b'^ — Ac is the discriminant. 

For example, 7 is a root of x^ — re — 1, which has discriminant D = 5. Therefore, 0(7) = l/v5- In §5, 
we study robustness of invariant circles for the real quadratic fields with the smallest discriminants, 
and correspondingly, the largest Diophantine coefficients. 

The Farey (or Stern-Brocot) tree is a recursive algorithm that can be used to generate the 
rationals between a given initial pair. These rationals can in turn be used to generate irrationals in 
a given algebraic ring. For convenience we write the rational number — as the vector {m,n), and 
will always assume that the greatest common divisor of the components is 1 . 

The root, or level zero, of the tree is a pair of rationals (mi, ni) and {ni2,n2) that are neighbors, 
i.e., 

min2 — Tiim2 = ±1. (34) 

We typically begin with the neighbors (0, 1) and (1,0), though sometimes we use other neighbors 
to zoom-in to some interval. The next level is obtained by the mediant operation © defined as 

(mi, ni) © (m2, n2) = (mi + m2, m + 712). (35) 

So, level one of the tree consists of the single vector (1, 1) = (0, 1) © (1,0), the daughter of the 
parent vectors (0, 1) and (1, 0). Note that the daughter is a neighbor to each of its parents. 

The {£ + 1)*** level is constructed by applying the mediant operation to each level i vector and 
its two parents; thus there are 2 rationals on each level. A binary tree is obtained by connecting 
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each vector at level ^ > 1 to its two daughters at level i + 1, see Fig. 18. There is a unique, finite 
path, the Farey path, from level one to every rational. This path is denoted by a sequence of Vs and 
r's denoting the left and right turns, respectively. Thus the path of (1, 1) is empty and, for example, 
p{5, 3) = rlr. Every irrational between the root rationals is uniquely the limit of an infinite path 
on the tree; for example, p{'j, 1) = rlrlrl . . . = {rl)°° , and p(l,7) = p(7~^, 1) = {lr)°° . According 
to a theorem of Lagrange, the Farey path for every quadratic irrational is eventually periodic and 
conversely every eventually periodic path converges to a quadratic irrational. 
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Figure 18: First five levels of the Farey tree construction. 

There is a simple relation between the continued fraction expansion of a number and its Farey 
path. Suppose that a positive rational has the continued fraction 

— = [ao,ai,... ,afc] = [ao,ai,... ,0^ - 1,1], 



where w.l.o.g. a^ > 1 for k >\. Then 






p(m, n)-s noiai^ao ,a,-l 



The continued fractions for the two daughters of — are obtained by incrementing the last entry by 
one in the two equivalent representations: [ao, ai, . . . , a^ + 1] and [oq, ai, . . . , afc — 1, 2]. 

Neighboring rational vectors can be used to generate irrational numbers in the ring 7j[uj]. To 
see this, note that (1,^;) is an integral basis for the ring, and every other integral basis will be of 
the form (?/i,7/2) = (l,w)M for some unimodular matrix M. Denote the rows of M by a pair of 
neighboring vectors (?7ii,ni), {m2,n2); then the irrational 

io = — = (36) 

r]2 ni + n2UJ 

is in Z[w]. The Farey path for uj' is essentially the concatenation of the paths of the rational 
and of UJ. More precisely, if {m2,n2) is a daughter of (77ii,ni), and its path ends in an /, then 
p{uj' , 1) = p{m2, n2)rp{l, uj) , while if it ends in an r, then p{uj' , 1) = p{m2,n2)lp{u}, 1). 

In §5, the 2+1 rational vectors up to level I on the Farey tree are used to construct 2 irrationals 
in the ring Z(a;) by applying the relation (36) to each pair of neighbors. 
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B Greene's Criterion 

Greene conjectured that periodic orbits in the neighborhood of an invariant circle should be stable. 
Indeed, a sequence of periodic orbits should limit upon a circle only if their Lyapunov exponents 
converge to zero. Conversely, if the limit of a family of periodic orbits has nonzero Lyapunov 
exponents, then the invariant circle should no longer exist. This is known as Greene's residue 
criterion. Although initially presented as a conjecture, some aspects of it have been proven [Mac92, 
FdlL92, DdlLOO]. 

A type-(7Ti,n) periodic orbit of / is an orbit of the lift, {xt+i,yt+i) = F{xt,yt), such that 
{xn, Vn) = (^^o + fn-, Do)- When / is area-preserving, the product of the multipliers of any periodic 
orbit must be 1. The stability of a periodic orbit is therefore completely determined by the trace 
of the Jacobian Df. Greene characterized this stability through a quantity he called the residue 

i?=i(2-trpr(xo,yo)). 

If < i? < 1, the orbit is linearly stable. 

Greene conjectured that an invariant circle with rotation number cv will exist if and only if the 
residues Rj of a sequence of type {ruj, Uj) periodic orbits that approximate the circle, 

lim — - = UJ, 

j-^oo Uj 

are bounded, 

lim sup \Rj\ < oo. 

j-s>oo 

Rational sequences that approximate to can be obtained by truncating either the Farey tree or 
continued fraction expansion of uj, recall App. A. To implement the criterion, we compute each 
periodic orbit using a root-finding method and continuation in e. This computation is simplified 
when the map is reversible, see App. C, since the computation can be reduced to a one-dimensional 
secant method. A threshold residue, R^'^, is selected, and the parameter value e*^, such that 
\Rj\ = i?*'* is found, also using the secant method. The choice of R^^ is arbitrary; however Greene 
found that for the golden mean, convergence was fastest with the choice R = 0.25 [Gre79]. This 
is the value we use. 

The approximation of Ecr can be improved by extrapolation. The renormalization theory implies 
that for noble circles of twist maps, 

for some constants A and 6 [Mac93]. This leads to a three-point extrapolation scheme: given three 
sequential approximants, we can compute the parameters 6, A, and Ecr- Extrapolation significantly 
improves the accuracy of the estimation of the critical parameter. For orbits up to period 30000, the 
values of e^J^ generally converged to 3 — 4 digits, but after extrapolation, the last several estimates 
of Ecr agreed to 4 — 8 digits. 

C Symmetries and Reversors 

Recall that a homeomorphism 5 is a symmetry of / if foS = Sof, and is a reversor if foS = Sof^^. 
The collection of reversors and symmetries of / is its reversing-symmetry group Q [LR98] . 
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An orbit is symmetric under 5 if it is invariant under S. One special class of reversors are orien- 
tation reversing involutions, and systems possessing such reversors have been called "R-reversible" 
[LR98]. Such a reversor has a one-dimensional fixed set Fix(5) = {z : S{z) = z}, and any sym- 
metric periodic orbit has a point on Fix(S) [Mac93]. Symmetric orbits can then be computed by 
a one- dimensional root-finding method along Fix(5) [Mei08]. It is for this reason that most of the 
studies of the breakup of invariant circles have been done for R-reversible maps and symmetric 
invariant circles. 

In this appendix, we recall some of the symmetries and reversors of the generalized standard 
map (1). Of course, any map / is a symmetry of itself. In addition, when / is a map on the 
cylinder, it commutes with integral rotations. Specifically, if 

Tm,nix,y) = ix + m,y + n), 171,11 eZ, 

then Tifi is a symmetry of /, or more properly of its lift, F, to the universal cover obtained by 
letting X G M. Since the composition of two symmetries is also a symmetry, F o Tmfi is also a 
symmetry of F. 

Certain special cases of (1) have additional symmetries and reversors. We say that a function 
g is odd about a if 

g{a + x) = —g{a — x), (odd about a). 

If the force in (1) is odd, then this map is reversed by the involution 

Si{x,y) = {2a-x,y + eg{x)). (37) 

Note that when g is periodic and odd about a, then it is also odd about ^ + a; the corresponding 
reversor is Ti^o ° Si. Every rotational invariant circle intersects Fix(5i) = {{a,y) : y G M}. The 
consequences for the conjugacy (7) are discussed in App. D. Cases studied in this paper with this 
reversor are shown in Table 5. 

Similarly, if the frequency map is odd about a point /3, then (1) is reversible under the 
involution 

S2{x,y) = {x,2P-y-eg{x)). (38) 

This reversor maps an invariant circle with rotation number w onto one with rotation number —oj. 
A consequence is that these circles have the same critical parameter set: ecr{—^) = £ct{'^)- 

Chirikov's standard map, (1) using (3) and (2), has both sets of reversors since il and g are 
both odd about zero (a = /3 = 0): it is "doubly reversible." Another way of viewing this is to note 
that this map has the inversion Ii{x, y) = (2a — x, 2/3 — y) as a symmetry. This is not independent 
of the reversors since I = Si o 82- 

An additional symmetry of maps, like Chirikov's map, with ^{y) = y is a translation symmetry 
in the momentum direction: / commutes with the vertical translation Tq^i. A consequence is that 
£cr{'-^) = £cri^ + n) for integer n. Combining this with the inversion implies that for Chirikov's 
map, £cr{^) = £cr(l " w); thus, one can limit the rotation numbers studied to cj G [0, ^]. The 
known symmetries of Chirikov's map are summarized in the first line of Table 5. 

The map (1) also has a reversor when both g and Jl are even, e.g., 

g{a + x) = g{a — x), (even about a). 
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When g is even about a and Q is even about /3, (1) has the reversor 

S3{x,y) = {2a-x,2P -y-eg{x)). 



(39) 



Note that this reversor is orientation preserving, and its fixed set is a point. This reversor maps 
an invariant circle with rotation number uj and positive twist onto a circle with the same rotation 
number but negative twist. Examples that we study, recall Table 5, include the nontwist maps 
with 0.2 ^ which is even about 0, and forces gi or g^, which are even about j and 



Finally, when Q is even and g has the odd-translation symmetry, g{x + 2] 
(1) commutes with the symmetry 

h{x,y) = {x + \,-y). 



4- 



-g{x), the map 



(40) 



This is the case for the standard nontwist map, i.e., using the Chirikov's force (2), and the frequency 
map (5). This symmetry was exploited in many studies of the breakup of shearless tori [SA98]. 

When none of the above symmetries of g and/or i} hold, then the map (1) is, as far as we know, 
not reversible, and its complete symmetry group is (/, Ti^o)- One example of this is due to Rannou 
[Ran 74]; however, while it is conjectured that this map is not reversible [Mac93, RQ92], as far as 
we know this question is still open. 

Map g 

Force Odd Even I il Odd Even I Reversors Symmetries 
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Table 5: Known generators of the reversing symmetry groups for various maps of the standard 
form (1), with forces 51-5, (2), (24), (25), (27), and (30), and frequency maps ili-s, (3), (5), and 
(26). The columns labeled "odd" and "even" give the points about which the functions are odd or 
even, if any. The "trivial" symmetries / and Ti^o are omitted. 



D Symmetries of the Conjugacy 

Symmetries or reversors act on the conjugacy k, (7), to produce additional conjugacies [AdlLPOS, 
OP08]. 
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Lemma 3. Suppose that k : S ^ M is a conjugacy (7) for f : M ^ M for a circle with rotation 
number oo. Then if S is a symmetry of f, k = S o k also solves (7), and if S is a reversor, 
k = S o k o R solves (7), where R{0) = —9. 

Proof. When S" is a symmetry, (7) implies that 

/oS'oA; = S'o/oA; = S'ofcoT(^. 

Thus k = S o k is also a solution to (7). When S is a reversor, then the inverse of (7) implies 

foSokoR = Sof^ okoR = Soko T_^ o R = S o k o R o T^. 

Thus k = S ok o R solves (7). D 

A symmetry or a reversor may map one invariant circle onto another, but if it maps a circle 
onto itself, then Lem. 3 and Lem. 1 imply that the conjugacy itself is symmetric up to a shift. One 
simple reversible example corresponds to the map (1) when g is odd, see App. C. 

Corollary 4. When (1) is a twist map, g is odd about a, k is a continuous conjugacy, and uj is 
irrational, then there is a (p such that kx — a is odd about tp, and ky is even about cp + i^co. 

Proof. For this case we use the reversor (37). The new conjugacy given by Lem. 3, 

k{e) = 5i o k{-0) = (2a - kx{-e), ky{-e) + sg{kx{-9)) , (41) 

is continuous and has degree-one, so that /c(S) is an invariant circle of / with rotation number 
UJ. When (1) is a twist map, it has at most one rotational invariant circle for each ui, therefore 
k(S) = k{S), and Lem. 1 thus implies that k{9) = k{9 + 2(p) for some (p. 
Consequently, the x-coordinate of (41) gives 

kx{9 + 2(p) = 2a-kx{-9). 

Setting S, = 9 + ip, then gives kx{p + S,) ~ ci = —{kx[p — C) — o/)', therefore, k^ — a \s odd about (p. 
Since the composition of a symmetry and a reversor is a reversor, a second reversor for (1) is 
f{Si{x, y)) = (2a — x + y, y). Composing (41), with / and using (7) this gives 

k{9 + u) = f{k{9)) = {foSi){k{-9)). 

Again using k{0) = k{9 + 2p), the y-component of the above equation gives ky{0+uj + 2ip) = ky{—9), 
which implies the evenness assertion. D 

We will use this result as a numerical check on our solutions of (7). While is seems difficult to 
check that a function is "odd" or "even" about some unknown point p, the Fourier coefficients of 
such a conjugacy, (8), must obey 

kxO = a- P, 
^{kxjc^'^n = 0, (42) 

9(^j,j-e2^"^'('^+5")) = 0, 
i.e., ip is determined by the average and the phases of the Fourier coefficients are related. 
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